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Abstract 

A Monte Carlo computer program (POEM) has been 
developed for the computation of photocurrents and 
photoemission in gamma (x-ray) -irradiated materials. 
The program has been used for application to the com- 
putation of radiation -induced surface currents on space 
vehicles and the computation of radiation-induced space- 
charge environments (IEMP) within space vehicles. The 
program has been successfully correlated with published 
expirmental data. 

POEM is a composite analytical -Monte Carlo com- 
puter program which calculates the photoemission from 
a gamma (x-ray) -irradiated slab of material. The pro- 
gram uses analytical routines to calculate photon flux 
attenuation, energy deposition and electron production 
within the slab, and uses a Monte Carlo routine to cal- 
culate the electron transport through the slab and the 
emission from the slab. The objective of using the 
analytical photon transport routine is to avoid the com- 
putation time-consuming complexity of a pure Monte 
Carlo model. When using a Monte Carlo photon trans- 
port routine, in order to produce a sufficient number of 
Compton scattering and photoelectric interactions to 
afford reasonable statistics for the initial energy-angle 
distributions of the electrons requires a great number 
of photon histories (typically several hundred thousand 
histories). To reduce the computation time required, 
an analytical routine has been incorporated in POEM to 
compute photon transport. 

The POEM program calculates the net photocurrent 
within a material and the electron emission current 
from a material including the energy flux, energy spec- 
trum, and angular distribution for an arbitrary incident 
photon spectrum. The emission can be calculated in 
either the direction of the incident photon flux or the 
reverse direction. In either case the program assumes 
normal incidence of the photon flux on the emitting slab, 
although with minor modification the program could 
treat an arbitrary angle of incidence. It is assumed 
that the time -dependence of the electron emission cur- 
rent is the same as the time -dependence of the incident 
photon flux, i.e. the model is quasistatic. 

The emission calculation includes the Compton 
electrons, photoelectrons (K and L), and Auger elec- 
trons (L) resulting from the scattering and absorption 
of the primary x-ray flux. The emission calculation, 
at present, does not include second generation (knock- 
on) electrons nor the electrons produced by the scat- 
tering and absorption of fluorescent or bremsstrahlung 
radiation. A study is planned to investigate techniques 
for the inclusion of these effects, including a more rig- 
orous treatment of the Compton scatter radiation. 


I. Introduction 

The description of the current density and energy- 
angle distribution of gamma ray or x-ray induced 
electron emission is required in the study of transient 
radiation effects on electrical systems and in the de- 
velopment of fast radiation detectors. A rigorous cal- 
culation of the electron emission requires both a Monte 
Carlo photon transport calculation of the electron 
source distributed throughout the irradiated matter and 
a Monte Carlo electron transport calculation of the in- 
tensity and energy-angle distribution of the multiple 
scattered electron flux at the emission surface. Be- 
cause of the low quantum efficiency for the production 
of primary electrons (Compton, Auger, and photo- 
electrons) by photons, it has been experienced that a 
very large number of photon histories are required in 
the Monte Carlo photon transport calculation to produce 
a statistically acceptable description of the electron 
source. Several hundred thousand photon histories 
have been found necessary for a typical incident photon 
environment, requiring on the order of 40 minutes 
computation time (using the ONE-D (*) code on the GE 
635 computer). The need exists for a less costly com- 
putational tool. A purely analytical electron emission 
model 3), on the other hand, while providing an 
inexpensive computation of the electron emission flux, 
does not properly treat electron multiple scattering 
and so does not provide a valid description of the 
angular distribution of emission. 

The POEM (photoemission) computer program 
was developed to provide a computational tool for the 
calculation of photon radiation-induced electron emis- 
sion which would be relatively inexpensive to use but 
would include a proper treatment of electron multiple 
scattering. POEM is a composite analytical -Monte 
Carlo program, using analytical routines to calculate 
photon flux attenuation and primary electron produc- 
tion in the material and a Monte Carlo routine to 
calculate the electron transport to the surface of the 
material. Whereas a typical electron emission cal- 
culation required 45 minutes computation time (on 
the GE 635) with a pure Monte Carlo code, the same 
calculation requires only 3 minutes computation time 
with POEM, without apparent, significant sacrifice in 
accurancy. 

POEM calculates the electron emission from a 
gamma or x-ray irradiated slab of material. The pro- 
gram assumes normal incidence of the x-ray flux on the 
emitting slab, although with minor modification the pro- 
gram could treat an arbitrary angle of incidence. The 
emission current, energy spectrum, and angular dis- 
tribution are calculated both in the direction of the in- 
cident radiation and in the reverse direction^ (Hence, 
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the program may be used to calculate the net primary 
photocurrent across any plane within the irradiated ma- 
terial.) It is assumed that the time -dependence of the 
electron emission current is the same as the time- 
dependence of the incident photon flux, i.e. the model is 
quasistatic. 

The emission calculation includes the Compton elec- 
trons, photoelectrons (K and L), and Auger electrons (L) 
resulting from the scattering and absorption of the pri- 
mary x-ray flux and (in an approximate treatment) the 
Compton scattered x-ray flux. The emission calculation, 
at present, does not include second generation (knock-on) 
electrons nor the electrons produced by the scattering 
and absorption of fluorescent or bremsstrahlung radia- 
tion. The electron transport calculation is strictly valid 
only for electron energies greater than the K absorption 
edge of the emitting material. The code seems to pro- 
duce reasonable results, however, down to somewhat 
lower energies. 

POEM has been evaluated through the comparison of 
computational results to published experimental data for 
electron emission induced by CuK^ x-radiation and Co^° 
gamma radiation, and through comparison to superflash 
x-ray-induced electron emission data obtained by the 
author. The comparisons to the published data are pre- 
sented herein. 


The total linear attenuation coefficient [i is then given by 


where T is the photoelectric cross-section and pis the 
density. Values for the fudge factor f have been selected 
by comparisons of exponential attenuation calculations to 
Monte Carlo calculations. 

II. 2 Electron Production 

The irradiated slab is divided into incremental slabs 
of thicknesses A x such that HA X << I, i. e. so that the 
photon flux variation through the zone is calculated in 
each zone for which there is a significant probability of 
escape, i. e. in each zone for which the maximum path 
length of the electrons produced in the zone is larger 
than the normal distance to the emitting surface. 
(Remember that the maximum electron path length is 
considerably greater than the extrapolated electron range 
range because of multiple scattering. ) 

The number of Compton electrons produced (per 
cm 2 ) in a zone, location x i? by photons of energy Ey is 

O A 

— Ax 

N c x i> = 1 ( E y > Kj) (1-e P ) 


II. Program Description 


II. 1 Photon Transport 


o 

The radiation intensity I (photons/cm ) at the point 
x of initiation of electron emission is calculated assum- 
ing a simple exponential attenuation of the intensity. 


I (x, Ey) = I o (Ey) e* +lx 


( 1 ) 


where I Q (Ey) is the incident intensity, x is the distance 
from the irradiated surface, and|i is the total linear at- 
tenuation coefficient. The assumption of exponential 
attenuation, of course, implies the complete absorption 
of a photon in any interaction. This is strictly true only 
of the photoelectric interaction. In a Compton scatter 
interaction, the photon recedes from the point of col- 
lision with reduced energy and altered direction. This 
Compton scattered radiation is treated by an approxima- 
tion: rather than reducing the photon flux by decreasing 
the energy of each Compton scattered photon while 
holding the number of photons constant (which would 
considerably complicate the calculation), the number of 
photons is reduced in each Compton interaction holding 
the photon energy constant. Mathematically this is done 
by introducing a Compton attenuation cross-section o T 
which is the sum of the Compton absorption cross- 
section and a fraction f of the Compton scatter 
cross-section a g , 


where O = q* + a s is the total Compton mass attenuation 
coefficient and the photon intensity I (E , x.) is given by 
equation The initial energy -angle probability dis- 
tributions of the Compton electrons are given by Klein- 
Nishina. 


The number of photoelectrons produced in the zone 
x by photons of energy Eyis 

--Ax 

N p (Ey, x.) = I(E y , x.) (1 -e P ) 

where T is the photoelectric mass attenuation coefficient. 
It is assumed that for all photon energies greater than 
the K edge the ratio of L to K photoelectrons produced is 
equal to the ratio at just above the K edge, i. e. 


T < E K - 

“k T(E K + C) - T(E k -€) 


( 6 ) 


where Ej^ is the K edge and € is an infinitesimal energy 
increment. The angular distribution of the photoelec- 
trons is given by Fischer and Sauter (4), the former for 
low energy photelectrons (j3<<l), the latter for high 
energy electrons (/3~ 1). 


The number of Auger electrons (L only) produced in 
the zone x ^ by photons of energy E^ is 


cr T = a + fa 
a s 


( 2 ) 


K 


N *<V*1> ■ < l -v % 


» X • ) 


(7) 
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where is the number of K photoelectrons produced 
and f^ is the K fluorescence yield given by.(^) 


III. Comparison of Computational Results 
With Experiment Data 


4 
z 

4 4 

z + 30 (8) 

where Z is the atomic number. The angular distribution 
of the Auger electrons is assumed isotropic. 

The total electron initial energy -angle probability 
distribution for each zone of the emitter is calculated by 
summing for each photon energy the Compton, photo- 
electric and Auger distributions and then integrating 
over the photon energy spectrum. 

II. 3 Electron Transport 


The transport of the electrons to the emitting sur- 
face is calculated using the POET ( 6 ) Monte Carlo elec- 
tron transport program which has been incorporated into 
POEM. POET is a condensed history Monte Carlo pro- 
gram. In the condensed history Monte Carlo method, 
the results of multiple scattering theory are combined 
with Monte Carlo techniques to form a consistent ap- 
proach to electron transport calculations in the media of 
finite dimensions. Since the large majority of electron- 
atom interactions involve small angel deflections and 
small energy loss, the condensed history random walk 
approach is applicable to electron transport ( it is not 
applicable to photon transport due to the large scattering 
angles involved in the majority of collisions). Rather 
than calculate the trajectories of the electrons as they 
undergo each coulomb interaction, a random walk com- 
putation is performed in such a way that each step takes 
into account the combined effect of many collisions. The 
electron history is divided into logarithmically-spaced 
energy intervals. For each step of the random walk the 
path length increment traveled by the electron, corres- 
ponding to the energy loss increment, is computed using 
the Bethe-M oiler stopping power. ^ By random 
sampling from the angular distributions calculated using 
the Goudsmitt-Saunderson theory of multiple scatter, 
the direction of travel of the electron after each random 
walk step is obtained. 

For each electron history the random walk com- 
putation is initiated by random selection of an initial 
location within the electron production zone (from a 
uniform distribution within the zone) and a random se- 
lection of the initial electron energy and velocity direc- 
tion (from the computed energy-angle probability dis- 
tribution for the zone). The random walk computation is 
continued until either the electron escapes from the 
emitting surface or until the residual path length of the 
electron is smaller than the normal distance to the dis- 
tance to the emitting surface (so that in this continuous 
slowing down approximation the electron has zero pro- 
bability of escape). 


III. 1 Energy Spectrum 

E. P. Denisov, et. al of Leningrad State Uni- 
versity published in 1963 the results of an expirment in 
which the structures of the x-ray photoemission energy 
spectra from aluminum, titanium, chromium, and iron 
photocathodes were measured. The photocathodes con- 
sisted of a thin layer of the test material (thicker than 
the maximum exit depth of the electrons) sputtered on a 
planar glass substratum at the center of an evacuated, 
spherical collector. The energy structure of the photo- 
emission was determined by measuring the collected 
current as a function of the bias potential between the 
photocathode and the collector. The photocathodes were 
irradiated by CuK a x-rays (8 kev). The experimental 
results for the four photocathodes are compared to the 
POEM calculations of the energy structure shown in 
Figure 1. The experimental and calculational results are 
compared by adjusting the arbitrary normalizations such 
that the primary maxima are equal. (No attempt was 
made to compare the gross emission currents since 
these data were net published by Denisov, et. al.). The 
agreement of the structures of the calculated and the 
experimental spectra is good except for the low energy 
spike in the experimental spectrum, which Denisov, et.al. 
interpret to be the secondary electron contribution to the 
emission spectrum. The three peaks in the energy 
spectrum due to the K photoelectron, L photoelectron, 
and Auger electron contributions to the spectrum can be 
seen in both the calculated and measured spectra. 

III. 2 Emission Current, Angular Distribution 

Ebert and Lauzon ^ of the Lawerence Radiation 
Laboratory published in 1965 the results of an experiment 
in which they measured the quantum efficiency (electrons 
emitted per incident photon) and the angular distribution 
of the gamma-ray induced electron emission from car- 
bon, aluminum, copper, cadmium, and lead targets. 

Two sets of target thicknesses were used: one set of 
targets with thickness equal to the range of a 1. 25 MeV 
photoelectron, and another set with thickness of one- 
tenth the range. The targets were irradiated by a col- 
limated beam of 1.25 MeV gammas from a Co™ source. 
The emission current was measured by a Faraday cup; 
the (cumulative) angular distribution was measured by 
varying the solid angle subtended by the Faraday cup. 

The Ebert and Lauzon results are compared to the 
POEM calculations for the quantum efficiency in Table 
1 and for the (differential) angular distribution shown in 
Figure 2. The agreement in both quantum efficiency and 
angular distribution is excellent. The agreement in 
quantum efficiency is not remarkable; Ebert and Lauzon 
had comparable success using an analytical model to 
calculate quantum efficiency. What is significant is the 


788 



References 



FARADAY 
TARGET CUP 
TEST CONFIGURATION 

QUANTUM EFFICIENCY { ELECTRON/PHOTON) 


TARGET 

THICKNESS (CM) 

POEM 

EXPERIMENT 

CARBON 

2.03 X 10" 1 

9.68 X IO“ 3 

9.14 X I0- 3 


2.26 X I0' 2 

2.57 X I0" 3 

2.46 X 10" 3 

ALUMINUM 

2.27 X IO _l 

7.96 X I0 -3 

7.97 X I0‘ 3 


1 .27 X I0 -2 

1.80 X I0“ 3 

1.99 X I0" 3 

COPPER 

5.24 X I0" 2 

6.10 X I0' 3 

6. 13 X I0‘ 3 


5.1 2 X I0' 3 

2.09 X I0" 3 

2.04 X IQ - 3 


TABLE 1. COMPARISON OF POEM RESULTS WITH 
EXPERIMENTAL DATA OF EBERT AND LAUZON 

agreement in the angular distribution. Ebert and Lauzon 
found that calculations of the angular distribution made 
using an analytical model, which assumed rectilinear 
electron transport, were in complete disagreement with 
the experimental data. Here, then, is the key advantage 
of the Monte Carlo calculation: analytical models which 
cannot properly treat electron multiple scattering, but 
rather assume rectilinear transport, provide no valid 
information as to the angular distribution of photoemis- 
sion. 

Conclusions 


Comparisons of computational results with experi- 
mental data have shown the POEM code to be quite 
adequate for the calculation of photon radiation-induced 
electron emission, at least in the limited ranges of 
photon energies for which test data were available for 
comparison. More test data, especially in the x-ray 
energy range, 1 keV to 1 MeV, are required for ad- 
ditional code verification. One shortcoming of the POEM 
code is its neglect of the secondary (knock-on) electrons. 
A study is planned of computational techniques for the 
inclusion of the secondary electrons in the emission 
calculation. 
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Figure la. Electron Emission Energy Spectrum - 

Aluminum Comparison of Calculations with 

Experimental Data (CuK O' x-Radiation 
Incident on Target) 
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Figure lb. Electron Emission Energy 
Spectrum - Titanium 
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Figure lc. Electron Emission Energy 
Spectrum - Chromium 
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Figure 2a. Electron Emission Angular Distribution - 
Carbon Target (Electrons Per Unit Solid Angle 
Per Incident Photon) (1.25 Mev Co^® Gamma 
Radiation on Target) 


Figure 2c. Electron Emission Angular 
Distribution - Cooper Target 



Figure 2b. Electron Emission Angular 
Distribution - Carbon Target 
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